%% Testing

% true label
Label_testing;
% estimated label
Label_est = [];
Label_est_b = [];

alpha = Data_testing*Weights;
Class_conditional = Sigma(alpha);
Label_est = round(Class_conditional);

% Switched to direct matrix operation by altering Sigma()
% for i = 1:N-2000
%     label_est_b = Sigma(alpha(i));
%     Label_est_b = [Label_est_b round(label_est_b)];
% end

Label_est_b = Label_est_b';

% Overlap = Label_est.*Label_testing;
% Correct = sum(Overlap);
dif = abs(Label_est - Label_testing);
Correct = max(size(find(dif == 0)));

Accuracy = Correct/(N-2000)
